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We review a novel approach to the superconductive proximity effect in disordered normal-super- 
conducting (N-S) structures. The method is based on the multicharge Keldysh action and is suit- 
able for the treatment of interaction and fluctuation effects. As an application of the formalism, 
we study the subgap conductance and noise in two-dimensional N-S systems in the presence of the 
electron-electron interaction in the Cooper channel. It is shown that singular nature of the inter- 
action correction at large scales leads to a nonmonotonuos temperature, voltage and magnetic field 
dependence of the Andreev conductance. 
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I. INTRODUCTION 



A superconductor in contact with a normal metal in- 
duces Cooper correlations between electrons in the nor- 
mal region, the phenomenon known as the proximity 
effect. Its microscopic origin lies in Andreev reflec- 
tion Q of an electron into a hole at the normal-metal- 
superconducting (N-S) interface. The probability of An- 
dreev reflection and thus the strength of the proximity 
effect is determined by the transparency of the N-S in- 
terface and the nature of electron propagation in the N 
part of the structure. Disorder in the normal conduc- 
tor near the N-S contact was shown theoretically |^-||] 
to increase considerably the effective probability of An- 
dreev reflection (see Ref. Q for a recent review from the 
experimental viewpoint). 

The standard semiclassical theory of N-S conductiv- 
ity based either on the traditional nonequilibrium 
superconductivity approach jl]] or on the scattering for- 
malism H[| , usually neglects interaction effects in the N 
part of the structure. However, in low-dimensional struc- 
tures, Coulomb interaction in the normal diffusive region 
gets enhanced Q , which may affect strongly the Andreev 
conductance and noise. 

In this paper we address the effect of interaction be- 
tween electrons in the normal part of an N-S structure 
on the charge transport through the system. To study 
a system with interaction a novel theoretical method 
should be developed since neither of the above-mentioned 
approaches can handle interaction corrections. Indeed, 
the scattering matrix formalism relying on the linear re- 
lation between the outgoing and incoming states is a 
priori a one-particle description. On the other hand, 
Larkin-Ovchinnikov kinetic equation jj] can be general- 
ized to allow for (at least some part of) interaction cor- 
rections, but its practical solution seems hardly possible 
beyond the first order of perturbation theory in interac- 
tion strength JlO 11|. 

An appropriate formalism had been developed in 
Rcf. p2| in the framework of the Keldysh action for dis- 



ordered superonductors |l3| . We start from the fully mi- 
croscopic Lagrangian describing interacting electrons in 
the diffusive conductor. Then, successively integrating 
over electronic degrees of freedom in the normal conduc- 
tor we end up with the Proximity Action, S prox [Qs, Qn], 
which is a functional of two matrices, Qs and Qn, de- 
scribing the states of the superconductive and external 
normal terminals of the N-S structure (cf. Fig. Q). Once 
the form of the Proximity Action is known, one can easily 
calculate the conductivity of the system, current noise, 
and, in principle, higher correlators of current and even 
the full statistics of transmitted charge |l^-|l5[ ■ 

The Proximity Action approach bears an obvious anal- 
ogy with the scattering matrix approach as both 
describe transport properties in terms of the character- 
istics of the terminals (stationary-state Green functions 
of the terminals Qs,n in the former versus asymptotic 
scattering states in the latter approach). In this respect, 
the Proximity Action method also shares the logic of 
Nazarov's circuit theory of Andreev conductance |I| . On 
the other hand, the Keldysh action approach is a natural 
generalization of the kinetic equation for dirty supercon- 
ductors for the case of fluctuating fields. The Larkin- 
Ovchinnikov kinetic equation then emerges as a saddle 
point equation for the Keldysh action flq] . 



V 




FIG. 1. A small superconductive island (S) of size 2d con- 
nected to a reservoir (R) through a tunnel barrier (I) and a 
dirty normal film (N) of size 2L 3> 2d. 
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As an application of the formalism, we will study 
charge transport in two-dimensional (2D) N-I-S struc- 
tures shown in Fig. |l| at low (compared to the S gap 
A) temperature and voltages, and arbitrary ratio t = 
Rd/Rt, where Rd and Rt are the resistances of the dif- 
fusive normal conductor, and of the tunnel barrier in the 
normal state, correspondingly. We will calculate the An- 
dreev conductance and noise of such systems in the pres- 
ence of Cooper interaction in the normal conductor mod- 
ified by the Coulomb interaction |16-18|, as a function of 
the "decoherence time" of an electron and the Andreev- 
reflected hole, h/£l*, where O* = max(T, eV, eDH/c). 
We will show that the Cooper interaction effects in 2D are 
substantially different from those in the ID case consid- 
ered in Ref. fTj}| , In particular, the lowest-order relative 
correction SGa/Ga scales as — X\n(L/d) and grows with 
the size of the system L, provided that £1* <C i?Th = 
HD/L 2 . We will sum up the main logarithmic terms 
of the order of [Aln(L / d)] n and find that Cooper repul- 
sion may lead to a nonmonotonuos dependence of G a on 
Rd/Rt and on the decoherence energy scale Jl*. 

In this paper we will not consider the weak localiza- 
tion |jl5^,|o) and interaction-induced corrections |9| to the 
sheet conductance, c, of the normal metal assuming that 
it is relatively large, g = (h/e 2 )a ^> \n(L/d). We will 
also not take into account the effect of the Coulomb zero- 
bias anomaly (ZBA) [ppl|] on the Andreev conductance. 
This effect was considered in |ll| within the lowest-order 
tunneling approximation. It was shown, similar to the 
Coulomb ZBA problem for normal conductors (cf., e. g., 
Ref. [Q), that it can be strongly suppressed if the bare 
long-range Coulomb potential is screened, e. g., by plac- 
ing a nearby metal gate. This is the situation we assume 
in the present paper. 

Another limitation of the present discussion is that we 
will consider the case of a 2D geometry of the current 
flow between the superconductive and normal electrodes 
of the structure, as shown in Fig. This will make 
possible to construct a unified functional renormalization 
group treatment that takes into account modifications of 
the proximity effect strength both due to multiple An- 
dreev reflections and due to Cooper-channel repulsion. 
Note that the sample geometry should not necessarily be 
symmetric as shown in Fig. ^: the only important feature 
of the geometry chosen is the 2D spreading of the current 
flow. 



II. MULTICHARGE PROXIMITY ACTION 

Here we present a brief overview of the Proximity Ac- 
tion approach to charge transport in N-I-S systems elab- 
orated in Ref. |L2]]. Eliminating all dynamic degrees of 
freedom in the normal diffusive conductor, one reduces 
the total action of the system to the form S prox [Qs, Qn], 
with the matrices Qs and Qat having the meaning of the 



steady-state Keldysh Green function of the S island and 
external normal terminal labeled by R in Fig. |], respec- 
tively. In the low frequency limit, ui <C E^h — D/L 2 , 
when diffusive electron motion spreads over the whole 
conducting region during one circle of the external volt- 
age, the resulting Proximity Action has the form 



SproxiQs, Qn] = -iit 2 g ^2 ^{QsQnY 

n=l 



(1) 



The action is known once the infinite set of parameters 
("charges") j n is specified. 

Remarkably, for noninteracting systems, knowledge of 
{"fn} is equivalent to knowledge of the whole distribution 
function V(T) of transmission coefficients of the system. 
To formulate this analogy, it is convenient to make a 
Fourier transform from {7,1} to a 27r-periodic function 
u(x) of an auxiliary continuous variable x: 



u(x) 



(2) 



The relation between the function u(x) and the gener- 
ating function of transmission coefficients introduced by 
Nazarov |^| (cf. also Ref. |23|) was established in Ref. |l2] ] 
by direct comparison: 



Ang u(x) 



T sin a; 



1 - Tsin^(o:/2) 



V(T) dT. 



(3) 



For interacting systems and systems with fluctuations 
p4j , Eq. (j|) may be interpreted as a generalization of 
the notion of transmission coefficients. 

Determination of {7™} is a separate task related with 
elimination of fast modes in the normal conductor. We 
will discussed it later, while now we will explain how to 
obtain transport properties from the Proximity Action 
(0). In order to study conducting properties of the sys- 
tem, we apply a bias voltage V to the superconducting 
and normal terminals of the structure and study its cur- 
rent response. The resulting expressions for the Andreev 
conductance (in units of e 2 /fi) and noise power take the 
form 111 



G A = 47rgu x (f ), 



(4) 



and 



(J W I- U ) = ^{(3-P S )*( W ) 

+ lp s [^(u;-2eV)+^(cj + 2eV)}}, (5) 

where fy(uj) = u> coth(o;/2T), and the superconductive 
noise function P5 is given by 



Ps = 1 



l^xxx ( 2 ) 

2«»(f ) ' 



(6) 
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The quantity F — (2/3)Pg gives the reduction of the shot 
noise between the normal and superconducting terminals, 
(I I) shot — F e{I), compared to its Poissonian value for a 
single electron tunneling, and can be identified with the 
Fano factor p5[ . 

For an N-S system, all quantities are expressed in terms 
of derivatives of u{x) at x — n/2. The counterparts of 
Eqs. (§-(§ for the case when the S island is in the nor- 
mal state can also be expressed in terms of the same 
function u(x), but with derivatives taken at x — 0. The 
conductance of the system G = Angu x (Q), while the 
current-current correlator is given by Eq. (^) with the 
replacements Ga —> G, 2eV — > eV, and P$ — > Pn = 
1 - 2ti ira:x (0)/'u x (0). 

In principle, one can go further and calculate the full 
charge transfer statistics in terms of the function u(x) 
following the derivation (ljjl^l in the scattering matrix 
technique. 



III. FUNCTIONAL RG 

We have seen that the set {~f n } or, equivalently, the 
function u{x) encodes all information about charge trans- 
mission through the system. Now we will discuss how to 
determine the charges 7^ for the system shown in Fig. [I]. 

The starting point is the microscopic Keldysh action 
for a S island in contact with a disordered metal, derived 
in Ref. |l3| . It can be represented as a sum of the bulk 
and boundary contributions. The bulk action, Sbuik, i s 
a functional of three fluctuating fields: the matter field 
Q(r,t,t') in the film, the electromagnetic potential, and 
the order parameter field used to decouple the quartic 
interaction vertex in the Cooper channel. Q(r,t,t') is a 
matrix in the direct 4 x 4-dimcnsional product K ® N 
of the Keldysh and Nambu-Gor'kov spaces, its average 
value giving the electron Green function in the Keldysh 
form, G(r, r). In the tunneling Hamiltonian approxima- 
tion, the boundary action, 



5, 



bound 



Stun [Qs, Q] 



iirGr 



TrQsQ(d), 



(7) 



describes an elementary tunneling process between the S 
island and the N metal. Here Gt — h/e 2 Rr is the di- 
mensionless tunneling conductance of the interface, and 
Q{d) is taken at the metal side of the N-S boundary. 

The next step is to eliminate degrees of freedom in 
the N film. 2D geometry of the system suggests that this 
procedure can be realized with the help of the Renormal- 
ization Group (RG) approach by successive integrating 
over fast modes in the bulk. As a result, the bound- 
ary action ([?]) which initially described single Andreev 
reflection process gets modified by multiple Andreev re- 
flections and acquires the form 



00 

■E 



7n(C)Tr[Q s Q(r)]' 



(8) 



where the logarithmic variable £ — 2 \n(r/d) — ln(u;d/fl), 
LOd = D/cP is the high-frequency cutoff, and £1 = D/r 2 
is the energy associated with the current scale r. At the 
energy scale aid, the multicharge action (|J) reduces to 
Eq. @j: 7 i(0) = a = G T /^g, and 7„> 2 (0) = 0, and, 
consequently, u(x, C = 0) = a sin a;. 

Under the action of the RG, the charges 7«,(C) are 
changing. In the zero-energy limit, when the Cooperon 
decoherence energy is smaller than the current RG en- 
ergy scale fl, both diffusons and Cooperons contribute to 
the RG equation which, written in terms of the function 
u{x,Q, reads [O 



+ uu x = — A(C) C) since. 



(9) 



The r.h.s. of the functional RG equation (f|) is due to 
the Cooper channel interaction with the scale-dependent 
constant A(£) given by |]l7| , ^3| 



A(C) 



Ad + A ff tanhA g C 
1 + — ^ tanhA s C 

A 9 



(10) 



where Ad is defined at the energy scale Lid- At ( > y/g, 
Coulomb repulsion in 2D drives A(£) toward the Fikcl- 
stein's fixed point \ g = XjtKyfg. In deriving Eq. (|^) we 
neglected the weak localization and interaction correc- 
tions to the film conductance g, making use of the small 
parameter Q/g <C 1. 

In the absence of Cooper channel interaction, Eq. (^) 
acquires the form of the Euler equation. In this case it 
is convenient to rewrite it in terms of the new variable 
t = Rd/Rt = oC an d the function U{x,t) = u(x,Q/a as 
Ut + UU X = 0, with the initial condition U(x, 0) = sinx. 
In this representation, all information about geometry of 
the system had gone, and the Euler equation for U(x, t) 
describes evolution of {7™} for any (not only 2D) system 
consisting of a tunnel barrier and a diffusive conductor 
in the zero-energy limit ]2^ , |26| , ^2] | . 

If the decoherence energy scale f2* > £/rh then Cooper- 
ons become inoperative at the latest stage of renormal- 
ization procedure, where £ > = ln(ajd/f2*). In this 
case the functional RG equation reduces to 



+ uu x = 0, 



(11) 



where u{x) = [u{x) — u{tt — x)]/2. The function u(x, £; C*) 
bears the knowledge of the decoherence scale through 
the initial condition u(x, C*;C*) = [ u ( x iC*) ~~ u ( 7r — 
x,(*)]/2, where u(x,(*) is the solution of Eq. (||) at the 
upper border of its applicability, £ = £*. Note, however, 
that crossover from the C < C* t° C > C* regimes is de- 
scribed by Eqs. (||), (]ll]) with the logarithmic accuracy 
only. 
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The RG equations (||) and (jll]) should be integrated 
until we reach the perfect boundary between the N film 
and the reservoir. Then, substituting Q(L) — Qn into 
the effective boundary action (@), we arrive at the Prox- 
imity Action (Q) with j n = 7„«Th), where Ctii = 
ln(uj d /E Th ) = 2ln(L/d). 

IV. ANDREEV CONDUCTANCE 

Now we proceed with the analysis of the functional RG 
equations derived in the previous section. We start from 
the zero-energy limit, f2* < Sxh, described by Eq. (||). 
In the absence of interaction one readily reproduces the 
known result 10,01 for the Andreev conductance: 



G 



(o) _ 



Gt 



sin0(t) 
l + isin9(t) 



= Gi 



tsinQ(t) 
1 + tsin0(t)' 



(12) 



where t = aCTh = Rd/Rti Rd = ln(L/e?)/27r<7, and the 
function 0(t) satisfies 0(f) = tcos0(i). 

It is often convenient to characterize the system with 
the help of the effective interface resistance Rt,cS defined 
formally through Ra = Rd + Rr,eS ■ For a noninteracting 
system, the ratio Rt^s/Rt as a function of t = Rd/Rt 
is shown in Fig. ^ by the dashed line. It behaves as t~ x 
for t <C 1 and saturates at 1 for t 3> 1. Correspondingly, 



G { a = G't/Gd for t <C 1, and G^' = G D for t > 1. 

At relatively small scales, the effect of Cooper interac- 
tion described by the r.h.s. of Eq. (|9|) can be considered 
perturbatively. The magnitude of the first-order inter- 



Co) 



action correction to the quasiclassical conductance (12), 
6Ga ~ — ACG? , grows 0] with the system size L and 

eventually becomes of the order of at £ ~ 1/A. Com- 
paring this scale to the scale £ ~ 1/a, corresponding 
to the crossover from the tunneling (t -C 1) to diffusive 
(i 3> 1) limits, one can distinguish between the regimes 
of strong, A»oe Gt/^9, and weak, A <C a, repulsion 
(still, in both cases A< 1). 




FIG. 2. Dependence of the effective interface resistance 
RT,es{i) / Rt vs. t = Rd/Rt for different values of A/a ob- 
tained by numerical solution of Eq. (^) for A(£) = const. The 
dashed line corresponds to the noninteracting case, A = 0. 



In general, the functional RG equation (^|) does not 
allow an analytical solution, and should be solved nu- 
merically. For clarity, we will assume hereafter that the 
interaction constant had reached the Finkelstein's fixed 
point X g = l/2n^/g, so that strong repulsion corresponds 
to a relatively weak tunneling conductance Gt <C 2-y/g, 
and vice versa. The effective interface resistance i?T,eff 
normalized to the tunneling resistance Rt as a function 
of t = Rd/Rt is plotted in Fig. || for different val- 
ues of the ratio X g /a. For the case of strong repulsion, 
X g ^> a, i?T,eff(i) very quickly (at t ~ a/X g <C 1) reaches 
its asymptotic value and saturates at RT, e s(t — oo) ~ 
(2X g /a)RT- The limiting value i?T,eff(°o) decreases with 
the decrease of X g /a up to X g /a ~ 1. At smaller X g /a, 
corresponding to the case of weak repulsion, Rt. c s(oo) 
starts to grow again, and reaches an asymptotic behav- 
ior Rt,cs(°o)/Rt ~ 1.19a/A g at X g /a <C 1. In this 
limit i?T,eff(i) reaches its asymptotic value at large scale 
t ~ a/Xg. 

The data of Fig. || demonstrate that in a 2D system 
the limits A — ► and Rd — > oo do noi commute: for any 
small but finite A 9 , i?T,eff(i)/ '-Rt will eventually (though, 
at very large i) deviate from the noninteracting depen- 
dence (dashed line in Fig. ||) and become large. 




fi*/E Th 

FIG. 3. Dependence of the Andreev conductance Ga(^*) 
(normalized to the zero-energy value) on the ratio fi*/-ETh 
for different values of t. For all plots, LUd/E^ = 400, A = A s , 
and Gt = 10 ^/g. 

At large enough temperature, voltage, or magnetic 
field, fi* = max(T, eV, eDH/c) ^> £rh, the resistance 
of the structure is given by p2] 



RA[t, t*) — Rd + RT,eS{t*), 



(13) 



where i* = u£* = {GT/^g) hi^d/fi*). Therefore, a 
nonmonotonous behavior of i?T, c ff(^*) will manifest itself 
in a nonmonotonous temperature, voltage and magnetic 
field dependence of the subgap conductance. The un- 
usual enhancement of conductivity with the increase of 
the decoherence energy scale f2* is most pronounced in 
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the limit of weak repulsion, Gt ^> ^\fg- Since the to- 
tal resistance is the sum of i?r,eff and Rd, cf. Eq. (13) 



the magnitude of the effect is determined by the ra- 
tio i?T,eff/ Rd which has a maximum at f x Gx/y/9- 
An example of such a nonmonotonous dependence of 
Ga(Q*) is shown in Fig. [|. The curves differ by the 
ratio t = Rd/Rt and correspond to Gt = 10-^/g (i. e., 
Xg/a = 0.2). 

We emphasize that the described nonmonotonous be- 
havior has nothing to do with usual "finite-bias" con- 
ductance anomaly in N-S structures with good inter- 
faces f2?],[ll]] , which occur at (T, eV ) « E Th even in the 
absense of interaction in the N region. In simple terms 
the origin of this new effect can be understood as follows: 
repulsion in the normal metal produces a superconduc- 
tive "gap function" in the normal conductor, Ajv, with 
the negative (compared to As in a superconductor) sign. 
Due to its opposite sign, An decreases the conductance 
of the structure, therefore any decoherence that reduces 
Ajv leads to the increase of the conductance. 




FIG. 4. Schematic map of the noise coefficient Ps as a 
function of the tunnel and diffusive resistances in the pres- 
ence of Cooper repulsion A = A 9 . Dark area refers to the 
tunnel limit Ps ~ 3, whereas light regions correspond to the 
diffusive regime with Ps ~ 1. 



V. NOISE 



Now we turn to the analysis of the noise function Pg. 
For a noninteracting system in the zero-energy limit it is 
given by |p8 12 



Pf (*) = ! + 



1 + G(t) tan 8(f) + 36(f) cot 6(f) 
2[l + 0(t) tan6(f)] 4 ' 



(14) 



Eq. ( J14| ) describes a crossover from the Poissonian (Pg = 
3) to the sub-Poissonian (Pg = I) |29| character of the 
noise as the system evolves from the tunnel (t <C 1) to 
the diffusive (t 3> 1) limits. 

In the presence of interaction, Pg can be estimated, 
qualitatively, by comparing the effective tunneling re- 
sistance, i?T,eff, with the diffusive resistance, Rd- For 
Rt.cH ^S> Rd, Ps ~ 3, while in the opposite case Pg w 1. 
If we again assume that A(£) = X g , then Pg will be- 
come a function of two parameters: t = Rd/Rt and 
Xg/a = 2^/g/GT- The boundaries between the regions 
with Pg = 3 and Pg = 1 on the plane (log t, log(C?T / \fg)) 
are sketched in Fig. [|. 

In the case of strong repulsion, Gt <C yfg, 
Rt,cS very quickly (at t ~ Gt/J~Q) saturates at 
RtM°°) = (^V9/Gt)Rt- Then one obtains |l| P s nt = 
Pn[Rd/Rt,ck(°o)], where P N (t) = 1 + 2/(1 + t) 3 is the 
noise function |j|,|l2| for the system with the N island. 
Thus, in the limit Gt y/g, the crossover between 
the tunnel and diffusive character of noise is shifted to 
t ~ ^g/G T » 1, see Fig. § 



In the limit of weak repulsion, Gt S> yfg, the situa- 
tion is more interesting. For t ~ 1, interaction corrections 
can be neglected and Pg is given by the noninteracting 
expression (|i~4|). So, at t ~ 1, Pg decreases from 3 to 
1, the corresponding boundary being shown in Fig. 
Later, at t ~ Gx/y/g 3> 1 (when Rd ~ h/e 2 y/g) inter- 
action corrections become relevant. In this region, i?T,eff 
is of the order of Rd, and one may anticipate that Pg 
will deviate from 1. For even larger t when resistance 
is dominated by the diffusive conductor, Pg will eventu- 
ally reduce down to 1. This crossover region is marked 
in Fig. H by the dashed lines. The function Pg in the 
crossover region t ~ Gt/ \fg 3> 1 obtained by numerical 
solution of Eq. ^ is plotted in Fig. | as a function of 
s = Xg( = (lnL/d)/TTy/g = 2(e 2 /K)Rd\/9- ft has a min- 
imum Ps = 0.99 at s = 0.40 and a maximum Pg = 1.28 
at s = 3.25. 




FIG. 5. Noise function Ps(s) vs. s = 2^/~§/Gd for the case 
of weak interaction, Gt 3> yf§- 
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Finally, we consider the fi* dependence of the noise 
power. In the case of weak interaction, Pg decreases 
(at t, S> 1) with the increase of the Cooperon deco- 
herence energy scale O*, as if the system is becoming 
more diffusive. This trend is opposite to what one has 
in the noninteracting case when the increase of 17* drives 
the system toward the tunnel limit, thus, increasing Pg. 
In the limit of strong repulsion, Gt *C ^fg, the zero- 
energy noise function Ps (t) exhibits a crossover from the 
tunnel to diffusive regimes at t ~ ^fg/Gr ^S> 1. Nev- 
ertheless, upon increase of 12* the function Pg(t,t*) re- 
mains t* independent down to much smaller values of 
t* ~ GtI \fg <§C 1 which correspond to the energy scales 
m(w d /12*) ^ 2ny/g. 



VI. CONCLUSION 

We developed the Proximity Action functional method 
able to describe quantum charge transport in mesoscopic 
superconductive structures in the presence of interaction 
in the normal metal. New method is applied to the study 
of the effect of repulsion in the Cooper channel upon An- 
dreev conductance and noise in 2D N-S structures. In- 
teraction corrections scale as <? -1 / 2 ln(L/d) and lead to 
nonmonotonous dependence of both the conductance and 
the Fano factor upon temperature, voltage and magnetic 
field. 
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